Eigenvector Approximation Leading to Exponential Speedup of Quantum Eigenvalue 
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We present an efficient method for preparing tfie initial state required by tfie eigenvalue approx- 
imation quantum algoritlim of Abrams and Lloyd 1]. Our method can be applied when solving 
continuous Hermitian eigenproblems, e.g., the Schrodinger equation, on a discrete grid. We start 
with a classically obtained eigenvector for a problem discretized on a coarse grid, and we efficiently 
construct, quantum mechanically, an approximation of the same eigenvector on a fine grid. We 
use this approximation as the initial state for the eigenvalue estimation algorithm, and show the 
relationship between its success probability and the size of the coarse grid. 

PACS numbers: 03.67.Lx, 02.60.-x 



Intuitively, quantum mechanical problems offer great 
potential for quantum computers to achieve large 
speedups over classical machines. An important prob- 
lem of this kind is approximation of an eigenvalue of a 
quantum mechanical evolution operator. In a recent pa- 
per 0, Abrams and Lloyd present a quantum algorithm 
for doing this. Their algorithm is exponentially faster 
than the best classical algorithm, but requires a good 
approximation of an eigenvector as input. In this pa- 
per we show how to obtain an approximation efficiently 
which is guaranteed to be good. 

The key component of the algorithm in 1] is quantum 
phase estimation, which is a method for approximating 
an eigenvalue of a unitary matrix 0]. We give a brief 
outline of this algorithm below. 

Let Q denote a 2™ x 2™ unitary matrix. We want 
to approximate a specific eigenvalue of Q. Phase esti- 
mation does this using the corresponding eigenvector as 
input. The algorithm in deals with the case when this 
eigenvector is not known exactly. In particular, consider 
a quantum computer consisting of three registers with a 
total of b + m + w qubits. The first b qubits are all ini- 
tially in the state |0). The second register with m qubits 
is initialized to some state {ip), which must approximate 
the eigenvector in question sufficiently well, as we will 
see. The last w qubits are work qubits for temporary 
storage. 

Since Q is unitary and therefore normal, the state {tp) 
can be expanded with respect to eigenvectors of Q. Omit- 
ting the work qubits, the initial state of the algorithm is 



|o)|V'>-|o)^rf„M, 



(1) 



where \u) are the eigenvectors of Q. Placing the first reg- 
ister in an equal superposition, using b Hadamard gates, 
transforms this state into 



(2) 



Next, powers of Q are applied to create the state 

2''-l 
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Since Q is unitary, its eigenvalues can be written as 
g27ri¥'„^ where (pu S R. We can assume that ipu S [O7 1) 
and consider the approximation of one of these phases 
instead of the approximation of one of the eigenvalues. 
Equation © is equal to 



^27ryip„ 
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It is easily seen that the inverse Fourier transform per- 
formed on the first register creates the state 



I 51 di^u, j)\j) 
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where 



I 1, 



2Vn ^ J 

2Vu = 3- 



(6) 

A measurement of the first register produces outcome j 
with probability 



J) 



and the second register will collapse to the state 
dug{(pu,j). > 



(7) 
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We remark that for the special case when the eigenvalues 
tpu can be represented exactly with 6-bits (i.e., 2^ipu is an 
integer), equation jSJ simplifies to 



(9) 
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Thus, when the eigenvalues are of this form an distinct 
a measurement of the first register will cause the second 
register to collapse exactly onto the corresponding eigen- 
vector. 

Recall that we are interested in approximating the 
phase that corresponds to an eigenvector \u'), that the 
state \ip) is an approximation of this eigenvector and 
that the eigenvalue is of the form e^'^*'''"' . For in- 
stance, one is often interested in the eigenvalue corre- 
sponding to the ground state. We define A{ipo,ipi) = 
mina;gz{|a; + ipi — ifo\}, ^o,Vi S R (i.e., the fractional 
part of the distance between i/jq and Lpi). Then a mea- 
surement of the first register produces an outcome from 
the set g ^ {j : A{j/2'>,ipu') < fc/2^ A: > 1} with 
probability 

> ^\du'g{'Pu',j)\'^ 
\du'\^ 



> \du'\- 
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(10) 



and when fc = 1 the probability that A(j/2^, (/?„/) < 2^^ 
is bounded from below by where the proofs of 

the probability bounds can be found in 0, 0] ■ Observe 
that must be chosen in a way that this probabil- 
ity is greater than i, which implies that |d„'| has to be 
sufficiently large. If we want to obtain an approxima- 
tion of (fiu' with accuracy 2~"' and probability at least 
Mii'P(l — e), equation l(Tn|l shows that this can achieved 
by choosing the number of qubits b in the first register 
to be 



log 1 



1 

Ye 



(11) 



The algorithm in [1| is based on the fact that quantum 
phase estimation can be used as an efficient subroutine to 
find eigenvalues. Consider a Hermitian operator of the 
form H = where each Hj acts only on a small 

number of qubits. Since H is Hermitian the operator 



Git) 



iHt 



is unitary and has the same eigenvectors as 



H. Using the technique in , G{t) can be approximated 
efficiently on a quantum computer and the approxima- 
tion is used as the unitary operator in the quantum phase 
estimation algorithm. (For more details regarding the 
efficient implementation of G see Q and the references 
therein.) 

The Hermitian eigenproblem described above is solved 
on a discrete grid. We are interested in the case when 
the grid is extremely fine. Clearly, a fine grid requires 
a large vector for the representation of the initial state 
of the algorithm. In general, it may not be possible to 
efficiently prepare an arbitrary quantum state in a space 
with a large number of qubits. However, in our case we 



will show a method for the efficient preparation of an 
initial state. 

Suppose we have an eigenvector for a coarse grid dis- 
cretization of the problem. We can assume that we ob- 
tained it classically because the size of the problem is 
small. Using this eigenvector, we efficiently construct 
an approximation to the corresponding eigenvector for 
a fine grid discretization of the problem. We use this 
approximation as the initial state of the eigenvalue ap- 
proximation algorithm. We describe our method for a 
one-dimensional continuous problem on the interval [0, 1] . 

Let H be an positive Hermitian operator, defined on 
a Hilbert space of smooth functions on [0,1]. Let Vk{-), 
k — 1,2,..., denote the eigenf unctions of H, ordered 
according to the magnitude of the corresponding eigen- 
values; and without loss of generality we assume that 



\vk{x)\'^dx = 1. 



(12) 



Suppose that Hn is a discretization of H with grid size 
hN ^l/{l + N). Let \ul^^), fc = 0,l,...,iV-l, denote 
the normalized eigenvectors of , ordered according to 
the magnitude of the corresponding eigenvalues. The ex- 
pansion of the fc-th eigenvector in the computational ba- 
sis can be written as 



N-l 



\J)- 



(13) 



Let \V^^'') = Yl^f^o^ Vk {{j + l)hN) \j) be the sampled 
version of Wfe(-) at the discretization points. Consider 
problems such that the eigenvector of interest satisfies 

IKIloo ^ supo<^<i \v'^^{x)\ ^ 0(1) and 



0{h%) 



(14) 



where q > is the order of convergence and |||^)|| 



N-l 



j=0 



for \X) — J2i=o ^ example, 



these conditions are satisfied when we are dealing with 
second order elliptic operators; see for the solution of 
eigenvalue problems. 



Now, assume that the eigenvector \U, 



of has 



been obtained classically 6]. This vector is placed in 
a logA^o qubit register. For N = 2'^No, we construct an 
approximation lujf^"^) of by appending s qubits, in 

the state |0), to \U^^°^) and then performing a Hadamard 
transformation on each one of these s qubits, i.e. 



|0) 



V2 



(15) 
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where f{j) = [j/2''\. The effect of / is to rephcate the The definition of \V^^^) and the fact that 



coordinates of luj.^"'') 2* times. We use |J7^"'') as in- 
put to the eigenvalue and eigenvector approximation al- 
gorithm. When the result of the algorithm is measured 
\uj.^^) will collapse onto a superposition of eigenvectors 
according to equation ((Sj). We show that the magnitude 
of the coefficient of \ujf^^) in this superposition can be 
made arbitrarily close to one by appropriately choosing 

Consider two different expansions of 1 1/^ '): 



N-l 

E 
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(16) 



(17) 



The first expansion is in the computational basis and the 
second is with respect to the eigenvectors of Hjy . We call 
Mfe^"*!^ probability of success. Equation H17(l can be 
rewritten as 



k I 



l^k 



(18) 

Taking norms on both sides and using ifT^ and lfTB|l gives 
the inequality 



N-l 

El 



(AT) 
k,k 



l=ik 



(Af)|2 



> 



E 

l^k 
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(19) 



We will now bound p9|) from above, and thus the prob- 
ability of failure. The definition of \tjjf^'') implies 



\U, 



{N 



Mifjj) + ^)hNo. 
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w, 



{N) 
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(20) 



where J2^_n^ lA^^^'^ 



= 0{h^j^) and Et^n' lAi, 



ANo) 12 _ 

2^0{h^^) by Applying the triangle inequality, we 

get 



It/, 




Vkiij + l)hN) 



VkiifU) + l)hNo) 



/2^|||y^(^°))|| 



1/2 



+ 0{h%). (21) 



(JV)\ 



imply that || |Vj, 
above is equal to 



felloe =0(1) 

N(l + 0{hN)). Hence, the sum 



N-l 



^E i^''((j' + i)M(i + o(/iw)) 



3=0 



-Vkiifij) + l)hr,o)il + Oih^M 



(22) 



Since Ufe(-) is continuous with a bounded first derivative, 
we have that 



(23) 



where xij = {j + l)hN and X2,j = (/(j) + l)/iWo, j = 
0, . . . , AT - 1. Clearly \x2,j - xi,j \ = 0{hNo)- Using 
(|23|l and the triangle inequality, we obtain from (|21|l that 



+0{h'' 



\u, 



\\\vr)\ 



N 



No' 



OihNo) 

(24) 



Hence, the probability of failure is bounded from above 
by 0{Nq '°'"{2,2(j}^^ depends only on the order of con- 
vergence to the continuous problem and the number of 
points in the classically solved small problem. We can 
select an Nq such that the probability of failure is less 
than 1/2, no matter how much larger N is. By choosing 
a large we can make the discretization error arbitrar- 
ily small. Equation H24|l implies that the probability of 

obtaining the eigenvalue e^'^"^'= with accuracy is at 
least A(i_o(jv- -'-{2.2,}))^ 

We remark that any classical numerical algorithm that 
computes an eigenvalue, satisfying a specific (nontrivial) 
property, oia NxN unitary matrix takes time ^l{N). For 
example, one may want to find the eigenvalue that corre- 
sponds to the ground state. This is true even if a matrix 
is sparse and regardless of whether the algorithm is de- 
terministic or randomized. It is merely a consequence 
of the fact that the algorithm needs to consider all the 
(nonzero) elements of the matrix and there are at least 
^}{N) of them. Alternatively, in the restricted case when 
the matrix is diagonal finding one of its elements is a 
problem at least as hard as searching an unordered list. 
The lower bound for searching yields the lower bound in 
our case. 

In conclusion, our method provides a highly efficient 
preparation of initial states for eigenvalue approxima- 
tion, requiring only a small number of Hadamard gates. 
Thus the algorithm of Abrams and Lloyd, using our ini- 
tial state, computes the eigenvalue exponentially faster 
than any classical algorithm. The method can be gener- 
alized to higher dimensional continuous problems. This 
will be the subject of a future paper. 
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